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FOREWORD 


This  report  was  prepared  by  the  Department  of  Mechanical  Engineering, 
Massachusetts  Institute  of  Technology,  under  Contract  No.  AF  l8(600)-957. 
The  contract  was  initiated  under  Project  No.  7351,  "Metallic  Materials," 
Task  No.  73521,  "Behavior  of  Jietals."  The  work  was  monitored  by  the 
Metals  and  Ceramics  Laboratory,  Directorate  of  Materials  and  Processes, 
Deputy  for  Technology,  Aeronautical  Systems  Division,  with  Mr.  D.  M, 
Forney,  Jr.  acting  as  project  engineer. 

This  report  covers  the  period  of  work  from  September  1959  to 
September  1961, 
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ABSTRACT 


nie  development  of  the  plastic  strain  in  a  V-grooved  flat  plate 
under  longitudinal  shear  was  followed  from  the  elastic  through  the  par¬ 
tially  plastic  to  the  fully  plastic  condition  for  a  non- strainhardening 
material.  Tie  region  of  plastic  flow  develops  monotonically.  Adjacent 
to  the  zone  of  deformation  in  the  fully  plastic  case  there  is  a  region 
where  limited  plastic  deformation  has  occurred. 

Tie  results  for  the  growth  of  the  plastic  zone  were  compared  with 
predictions  based  on  the  elastic-plastic  solution  for  an  infinite  plate 
and  the  elastic  solution  for  a  finite  plate.  Agreement  is  good  at  low 
stress  levels.  At  high  stress  levels,  a  relatively  simple  empiricsLl  e- 
quation,  satisfying  overall  equilibrium,  is  proposed.  Predictions  based 
on  elasticity  theory  alone  are  shown  to  be  seriously  in  error. 
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I.  INTRODUCTION 


For  ein  understanding  of  fatigue  it  is  desirable  to  know  the  dis¬ 
tribution  of  plastic  yielding  at  the  tip  of  a  groove  or  growing  fatigue 
crack.  Since  problems  in  longitudinal  shear  are  more  readily  solved 
than  others,  it  was  decided  to  study  a  circumferentially  notched,  thin- 
walled  tube  loaded  in  torsion,  which  could  be  idealized  as  a  grooved 
plate  loaded  in  longitudinail  shear,  as  shown  in  Fig.  1.  The  plate,  con¬ 
taining  an  infinite  V-groove  in  its  face,  is  loaded  with  a  shear  stress 
'T'co  parallel  to  the  groove.  Ihe  plate  thickness  is  t,  the  groove  depth 
is  c,  and  the  semi-groove  angle  is  Ct  .  The  material  is  elastic-plastic, 
non- strainhardening  with  a  yield  stress  in  shear  of  k. 

II.  ELASTIC-PLASTIC  ANALYSIS 


Consideration  of  the  case  of  an  infinitely  thick  plate,  given  by 
Hult  and  Mcdintock  (195^),  or  of  the  limiting  case  of  torsion  as  given 
by  Prager  and  Hodge  (l95l)>  leads  one  to  postulate  displacement  only  in 
the  direction  of  the  axis  of  the  groove.  Uie  equations  of  equilibrium 
are  satisfied  if  the  components  of  stress  are  set  equal  to  derivatives 
of  a  potential  function: 


d0 

dy 


a  <P 


(1) 


In  the  elastic  region  Hooke's  Law  and  the  dependence  of  strains  on  the 
one  component  of  displacement  then  lead  to  the  familar  Laplace  equation 
for  the  stress  potential 

subject  to  the  boundary  condition, 

^  (>« »y)  +  1^ 

integrated  around  the  boundary  of  the  elastic  region* 


0 


(2) 


Equation  2  with  its  associated  boundary  condition,  Eq.  3/  applies 
only  in  the  region  where  Hooke's  Lai;  is  satisfied  and  can  be  solved 
only  with  a  knowledge  of  the  position  of  the  elastic-plastic  boundary. 
(More  precisely,  it  is  not  the  plasticity  but  the  non-linearity  of  the 
stress  strain  relation  which  is  critical  in  determining  this  boundary. 
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but  the  two  conditions  are  identical  for  the  material  considered  here. ) 
The  need  to  know  the  elastic-plastic  boundary  is  avoided  by  transforming 
from  space  coordinates  to  stress  coordinates.  Since  the  Jacobian  of 
and  with  respect  to  x  and  y  is  non-zero,  the  continuity  and  equili¬ 
brium  equations  can  be  written  as 


The  second  of  these  will  be  identically  satisfied  if  the  coordinates 
are  set  equal  to  the  partial  derivatives  of  a  potential  with  respect 
to  the  components  of  stress: 


( 


a 


From  Eq. 
equation, 


k,  the  coordinate  potential  then  must  satisfy  Laplace's 


(5) 


^  Txz^  dTyz^  ^  * 


(6) 


subject  to  the  boundary  condition 

'Tnz.Tyt 


■/( 


dTrz 


xz 


(7) 


This  transformation  and  the  resultin'  rc.ndary  conditions  are  shown 
in  Fir.  :  .  In  interratinr  aro'und  the  elastic-plastic  bo'indary  it  must 
be  noted  that  the  resultant  stress  is  equal  ts  the  yield  stress  and  is 
normal  to  a  posit: cn  vector  from  the  tip  cl  tre  groove.  The  entire 
bo'undary  except  the  magnitude  of  the  stress  at  the  free  surface  op¬ 
posite  the  groove  (  9^  )  is  .uiovn.  rC'.'ever,  :‘rcr.  the  continuity  of 
stress  :n  space  it  can  be  shown  that  {d^/dTxz  )  is  continuous 

in  the  7rz  “  direction  across  the  toint  (  Tiz  »  Ot^YZ-  7b),  which 
is  s deficient  information  to  uniquely  determine  a  solution. 


Hue  to  this  indeterminate  bo’uncary  ccnciticn,  the  author  was  not 
able  to  find  a  closed  fom,  analvtic  solution  to  the  prcclem.  There¬ 
fore  the  scluPion  carried  cut  by  machine  computation^  using  relaxa- 


^  These  com pu tat ions  were  performed  In  part  at  the  K.  I.  T.  Compu¬ 
tation  Center,  Cambridre,  I'ess.,  and  ccr.pleted  at  the  University  Mathe¬ 
matical  Laborator:/,  Cambridge,  England. 
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tion  techniques  plus  a  provision  for  satisfying  continuity  of 
across  7yz  =  Tb  •  To  simplify  computation,  the  sector-shaped  area  of 
the  stress  plane  was  mapped  conformally  into  a  semi-infinite  strip#  Since 
the  zero  stress  point  of  the  stress  plane  is  a  singular  point  of  this 
transformation,  the  boundary  there  was  replaced  by  a  small  circle  of 
radius  )  and  ^  was  simply  called  zero  on  this  circle#  The 

number  of  cycles  to  obtain  proper  convergence  and  the  accuracy  of  the 
finite  difference  approximations  were  determined  as  outlined  by  Crandall 
(1956)#  The  degree  of  convergence  was  checked  by  performing  extra 
cycles  in  one  case#  The  slopes  of  the  potential  function  are  accurate 
to  about  5/^  and  are  on  the  low  side  of  the  correct  value.  The  results 
of  these  computations  are  shown  in  Figs.  3  through  8. 

Hult  and  McClintock  (1956)  have  sho\m  that  in  the  plastic  zone,  the 
strain  consists  of  a  single  shear  component  normal  to  the  position  vector, 
f#  This  shear  component  is  given  in  terms  of  the  yield  strain,  Yyk/Q, 
and  the  distance  to  the  plastic-elastic  boundary  (  R  ),  rneastured  through 
the  point  in  question,  by 


y»z  =*  /y  • 


(8) 


Iherefore  since  the  elastic-plastic  boundary  is  known,  the  strain  dis¬ 
tribution  in  the  plastic  zone  can  be  determined. 

It  can  be  seen  from  Fig.  6  that  for  any  element,  the  strains  are 
monotonically  increasing  functions  of  the  ap^plied  load.  Hierefore  the 
same  result  is  obtained  for  either  non-linearly  elastic  or  elastic- 
plastic  material. 

Once  the  fully  plastic  state  is  reached,  all  further  deformation 
takes  place  on  the  plane  y  =  0,  as  predicted  by  a  rigid-clastic  analysis. 
It  is  observed  that  adjacent  to  the  plane  of  deformation  in  the  fully 
plastic  state,  a  region  of  limited  deformation  exists  as  suggested  by 
Hill  (1950). 


III.  CO'IPARISON  V/ITH  ELASTIC  ANALYSIS 


Since  numerical  computations  are  time  consuming  and  do  not  afford 
a  cCJmpact  summary  of  results,  it  is  well  to  inquire  how  closely  these 
results  might  be  emticipated  from  an  elastic  analysis  coupled  with  the 
elastic-plastic  analysis  of  Hult  and  McClintock  (1956)  for  a  plate  of 
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infinite  thiciness. 


As  shorn  by  iAlllai’-s  (l^'57)  ^'or  zero-anj:;le  notches  under  tension  or 
transverse  shoar^  the  stress  distribution  near  the  tip  of  the  notch  in 
a  body  of  aroltro.ry  shap.e  can  be  characterized  by  a  sin::le  pararrieter. 
oinilarly,  for  :i  zero-an  :le  jroove,  or  crach,  ’under  lon^'^itudinal  shear, 
halt  and  Acflintoc*'  fo'und  the  elastic  stress  distribution  to  be 

=  Kjr-'’/^cosf  ;  Tp^  =  Kjr .  O) 

for  ih-c  *:nfinito  lute,  the  constant  .0:13  .^iven  In  teia^.s  of  the  stress 
at  iir.'nlty  (  7^  )  ■■nd  tl.e  crac..  de't?!  (c)  by 

K300  =■  750(0/2)'/^,  (10) 

for  a  :  ater’'.l  *.'1*::  y!^f  d  strenjth  'n  shear  of  I:,  the  radius  of  the 
plastic  uens  a  n  b*:  to  r.d  fro.-  the  s  se  of  a  crach  in  an  infinite  solid 
to  be  ,;iven  'r.  of  the  cja::t’c  :tr'"S3  intensity  factor,  of  the 

el''*stlc  sc"  '.  :'en  * ; ' 

R  •  2Ksyk*  =  C  (7i,/k)®  , 


n.e 


factor  for  t:,e  'late  of  finite  thick¬ 
et  Icn  s'lr.ilar  to  that  used  by  'lester- 


ot 


■-  si  a  :;cr:^a  of  cr" z  In  a  plrtc  .mder  tension. 
3  '  sor*'e St*  on  for  ifatc  thlciuieso  Identical  to 
'd,  '.'h-ch  .*s  '"Iso  reco.s  ended  by  the  A^i!l  (l9i'0) 


(12) 


exact,  in  the 
stress  exists 
analysis,  but 


ite  3*. 


cv,  re  the  a  ’'lied  streso  's  increased 
'  ‘  “f  th^n  ’:  -*en  by  hcs.  11  and  1?.  An 
r.  f.atad  by  noting  that  t'hen  the  r.ini- 
;.lc,  squilii  r ! vr.  requires  that  the  ap- 
r:'’.s  can  be  achieved  by  adding  a  cor- 


Choosing  the  exponent  on  the  stress  term  to  be  5^ 
give  the  proper  limit  gives 


(t-c)(1  - 


2(t-c) 


t<in  ^  ) 


and  the  constant  C  to 


(iM 


The  results  of  the  approximate  relation  given  by  Eq,  13  are  plotted  in 
Figs*  3  through  5  and  jive  satisfactory  agreement  with  the  more  exact 
numerical  solution.  It  is  also  of  interest  to  canpare  these  results 
with  the  assumption  that  the  plastic  2:one  extends  to  the  point  T/here  the 
stress  calculated  on  an  elastic  basis  reaches  the  yield  stress,  as  re- 
coi:u:ended  by  the  AS’H'l  (l9'^0).  In  rxahlnj  the  correction  for  the  finite 
width  of  a  plate,  they  sujjest  that  the  effective  crack  length  shoiiLd 
be  tal'ien  to  be  the  actual  craci:  length  ‘1*.  s  the  radius  from  the  tip  of 
the  crack,  to  the  elastic-plastic  bovn^ar;^’’  directly  ahead  of  the  crack, 
c  +  Ihese  ideas  lead  to  the  equation 

As  shown  in  Fijs.  3  tl'roujh  5^  tr.is  approximation  is  unsatisfactory.  In 
the  first  place,  at  low  stress  levels  it  was  3110^^  by  Ilult  and  IlcClintock 
(l95o)  that  this  type  of  aprroxir  umion  gave  res’flts  lov  by  a  *actor  of 
tv;o,  Cn  the  otherhand,  at  high  stress  levels  the  eqi::ation  is 
unsatisfactory*'  because  it  implies  that  wren  the  plastic  zone  reaches  all 
the  way  across  the  plate,  the  load  carrying  ca2)acity  has  dropped  to  zero. 
Even  ;n.thout  tl  e  appearance  of  the  term  for  the  radius  of  the  plastic 
zone  in  the  correction  for  finite  pla-e  'width,  the  equation  fails  to 
satisfy  equilibrium  in  the  fully  plasvic  case.  Since  this  objection 
holds  in  the  case  of  a  sheet  under  tension  well  as  a  plate  under  shear, 
the  tensile  analog  of  hq.  15  is  to  be  preferred. 

A  further  disadvantage  to  the  above  assuir.ption,  that  the  boundary 
of  the  plasmic  zone  is  the  point  at  v/hich  the  elastic  stress  distribu¬ 
tion  reaches  the  yield  stress,  is  seen  on  noting  that  such  an  assump¬ 
tion  would  indicate  plastic  flow  occurrirg]  along  the  fianl:  of  the 
groove,  \/hercas  as  shov/n  in  F-gs.  6  througl:  3,  the  plastic  zone  actually 
lies  entirely  ahead  of  the  craci:. 


5 


IV.  CONCLUSIONS 


1)  In  the  fully  plastic  case^  the  strain  in  a  singly-grooved 
flat  plate  under  longitudinal  shear  is  concentrated  along  the  minimiam 
section.  As  the  plastic  flow  develops,  however,  small  scale  plastic 
strain  occurs  in  a  monotonically  increasing  region  on  either  side  of 
the  plane  of  minimum  cross  section. 

2)  The  extent  of  the  plastic  zone  ahead  of  the  groove,  R^,  is 

given  in  terms  of  the  plate  thickness  t,  the  crack  length  or  groove 
depth  X,  the  nominal  applied  stress,  and  the  yield  stress  in 

shear,  k,  by  the  equation 


where  the  constant  C^  chosen  to  satisfy  equilibrium  in  the  fully 
plastic  case,  is 

If-)  . 

3)  A  co’“responding  equation,  derived  by  a  r.ethod  recommended  by 
the  ASTK,  is  in  error  by  a  factor  of  two  at  low  stress  levels,  and  is 
even  qualitatively  incorrect  at  high  stress  levels. 

4)  Within  the  plastic  zone,  the  strain  is  given  in  terms  of  the 
radius  to  the  elastic-plastic  boundary  by 
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FI6.3  EXTENT  OF  PLASTIC  ZONE  VS.  APPLIED  LOAD, 
WITH  APPROXIMATE  EQUATIONS  (c/t  =  1/4). 
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APPLIED  STRESS/ FULLY  PLASTIC  STRESS  t/k(t-c) 


FIG. 4  EXTENT  OF  PLASTIC  ZONE  VS.  APPLIED  LOAD, 
WITH  APPROXIMATE  EQUATIONS  (c/t  =  l/2) 
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EXTENT  OF  PLASTIC  ZONE/EXTENT  OF  FULLY  PLASTIC  ZONE.  IRq/I-c) 


APPLIED  STRESS/ FULLY  PLASTIC  STRESS  (root/k(t-c) 


FIG.  5  EXTENT  OF  PLASTIC  ZONE  VS.  APPLIED  LOAD, 
WITH  APPROXIMATE  EQUATIONS  (c/t  =  3/4) 
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